Chromoelectric flux tubes and coherence length in QCD 
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The transverse profile of the chromoelectric flux tubes in SU(2) and SU(3) pure gauge theories 
is analyzed by a simple variational ansatz using a strict analogy with ordinary superconductivity. 
Our method allows to extract the penetration length and the coherence length of the flux tube. 

PACS numbers: 11.15.Ha, 12.38.Aw 



I. INTRODUCTION 

The presence of chromoelectric flux tubes in QCD 
vacuum is a clear signal of color confinement [l|, 
Monte Carlo simulations of lattice QCD can produce a 
sample of vacuum configurations, thus allowing a thor- 
ough nonperturbative study of tube-like structures that 
emerge by analyzing the chromoelectric fields between 
static quarks [3HI9I]. A direct consequence of the tube- 
like structure of the chromoelectric fields between static 
quarks is the linear potential and hence the color con- 
finement. 

A striking physical analogy exists between the QCD 
vacuum and an electric superconductor. As conjectured 
long time ago by 't Hooft [2(| and Mandelstam plj . 
the vacuum of QCD could be modeled as a coherent 
state of color magnetic monopoles, what is well known 
as dual superconductor 22] . In the dual superconductor 
model of QCD vacuum the condensation of color mag- 
netic monopoles is analogous to the formation of Cooper 
pairs in the BCS theory of superconductivity. Even if the 
dynamical formation of color magnetic monopoles is not 
explained by the 't Hooft construction, there is a lot of 
lattice evidences [23l - l3l| for the color magnetic conden- 
sation in QCD vacuum. It should be recognized [32j that 
the color magnetic monopole condensation in the confine- 
ment mode of QCD could be a consequence rather than 
the origin of the mechanism of color confinement, that 
actually could be originated from additional dynamical 
causes. Notwithstanding the dual superconductivity pic- 
ture of the QCD vacuum remains at least a very useful 
phenomenological frame to interpret the vacuum dynam- 
ics. 
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In the usual electric superconductivity tube-like struc- 
tures arise 33[ as a solution of the Ginzburg-Landau 
equations. Similar solutions were found by Nielsen and 
Olesen 34] in the case of the Abelian Higgs model, where 
they showed that a vortex solution exists independently 
of the type I or type II superconductor behavior of the 
vacuum. In previous studies [T2TH6I l35j performed by 
some of the present authors, color flux tubes made up 
of chromoelectric field directed along the line joining a 
static quark- antiquark pair has been investigated, in the 
cases of SU(2) and SU(3). 

In the present work we would like to push forward the 
analogy with electric superconductivity and exploit some 
results [36[ in the superconductivity to further extract 
information from flux tube configurations in SU(2) and 
SU(3) vacuum. The method and the numerical results 
for both SU(2) and SU(3) are reported in Section [TTJ In 
Section IIIII we check the scaling of the penetration and 
coherence lengths for both SU(2) and SU(3), and com- 
pare with previous studies. In Section IIVI we critically 
discuss the contribution of the longitudinal chromoelec- 
tric field to the string tension. Finally, in Section IVl we 
summarize our results and present our conclusions. 



II. CHROMOELECTRIC FLUX TUBES ON THE 
LATTICE 

To explore on the lattice the field configurations pro- 
duced by a static quark- antiquark pair we exploit the 
following connected correlation function 0, H, H3, HH 
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where Up = U^ u (x) is the plaquette in the (p,v) plane, 
connected to the Wilson loop W by a Schwinger line L, 
N is the number of colors (see Fig. 1 in Rcfs. [HI, HE]). 
The correlation function defined in Eq. ([1} measures the 
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field strength. Indeed, in the naive continuum limit 
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where ( ) q q denotes the average in the presence of a 
static qq pair and ( )o is the vacuum average. According 
to Eq. ([2]), we define the color field strength tensor as 
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By varying the distance and the orientation of the plaque- 
tte Up with respect to the Wilson loop W, one can probe 
the color field distribution of the flux tube. In particular, 
the case of plaquette parallel to the Wilson loop corre- 
sponds to the component of the chromoelectric field lon- 
gitudinal to the axis defined by the static quarks. In pre- 
vious studies the formation of chromoelectric flux tubes 
was investigated in SU(2) lattice gauge theory [Tol fT3 - fl6j 
and in SU(3) lattice gauge theory [35| by exploiting the 
connected correlation function Eq. (JTJ . It was found that 
the flux tube is almost completely formed by the longitu- 
dinal chromoelectric field, E\ , which is constant along the 
flux axis and decreases rapidly in the transverse direction 
x t . By interpreting the formation of chromoelectric flux 
tubes as dual Meissner effect in the context of the dual 
superconductor model of confinement, the proposal was 
advanced [l(3, Il2l - fl6| to fit the transverse shape of the 
longitudinal chromoelectric field according to 

Ei(x t ) = -^fi 2 K (^x t ) , x t > . (4) 

Here, Kq is the modified Bcssel function of order zero, <f> is 
the external flux, and A = l//z is the London penetration 
length. Equation (TJ| is valid in the region x t ^>^,£, being 
the coherence length which measures the coherence of the 
magnetic monopole condensate (the dual version of the 
Cooper condensate). In fact, we expect that Eq. ((4]) gives 
an adequate description of the transverse structure of the 
flux tube if A 3> £. This means that Eq. Q should be 
valid for k 1 (type II superconductor), where k is the 
Ginzburg-Landau parameter, 
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However, several numerical studies [1, l39U48j in both 
SU(2) and SU(3) lattice gauge theories indicated that 
the vacuum behaves like an effective dual superconduc- 
tor which belongs to the borderline between a type I and 
type II superconductor with k ~ 1. Thus, we see that 
Eq. (fJJ is no longer adequate to account for the trans- 
verse structure of the longitudinal chromoelectric field. 
Remarkably, it turns out that we may re-analyze our lat- 
tice data for chromoelectric flux tubes by exploiting the 
results presented in Ref. 36J where, from the assumption 



of a simple variational model for the magnitude of the 
normalized order parameter of an isolated vortex, a sim- 
ple analytic expression is derived for the magnetic field 



and supercurrent density that solve Ampere's law and the 
Ginzburg-Landau equation. In particular, the transverse 
distribution of the magnetic filed reduces to the London 
model results outside the vortex core, but has the added 
advantage of yielding realistic values in the vortex core 
vicinity. Accordingly, from Eq. (4) of Ref. [361 ] we derive 

4> 1 K (R/X) 
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where is a variational core radius parameter found to 
be [3(| of the order of £. Equation © can be written as 
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By fitting Eq. ([5]) to our flux tubes data, we may obtain 
both the penetration length A and the ratio of the pen- 
etration length to the variational core radius parameter 
A/ft,. It is worth to recall that, by means of Eq. (JHJ, 
we can extend our fit up to xt — 0. Moreover by using 
Eq. (16) of Ref. (36|, we may also obtain the Ginzburg- 
Landau k parameter, 
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with K\ the modified Bessel function of order 1. The 
coherence length f is obtained from Eqs. (O and (fTU|) . 
Our data for chromoelectric fields between static quark- 
antiquark sources have been obtained through the con- 
nected correlation function Eq. (fTJ). In order to reduce 
the quantum fluctuations we adopted the controlled cool- 
ing algorithm. It is known [49( that by cooling in a 
smooth way equilibrium configurations, quantum fluctu- 
ations are reduced by a few order of magnitude, while 
the string tension survives and shows a plateau. We 
shall show below that the penetration length behaves in 
a similar way. The details of the cooling procedure are 
described in Ref. [16| for the case of SU(2). Here we 
adapted the procedure to the case of SU(3), by applying 
successively this algorithm to various SU(2) subgroups. 
The control parameter S was fixed at the value 0.0354, as 
in Ref. [l6[ . As described in Ref. [35] , in the construction 
of the lattice operator given in Eq. (fTJ) we have consid- 
ered also noninteger distances to check the restoration of 
the rotational symmetry on our lattices. 



A. SU(2) data 

We analyzed our lattice SU(2) data collected for three 
different values of f3, namely B = 2.52,2.55,2.6 (for fur- 
ther details we refer to Ref. [35]). Indeed, we find that 
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Eq. © is able to reproduce the transverse distribution of 
the longitudinal chromoelectric field in the whole region 
Xt > 0. An example of the of effectiveness of Eq. ((HJ to 
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FIG. 1. SU(2): Ei versus x t for /3 = 2.52 after 8 cooling steps. 

fit all data for the transverse distribution of the chromo- 
electric field down to Xt = is given in Fig. [T] where 
we also display the points calculated at nonintcger dis- 
tances, which were not included in the fit. We see that 
there are slight deviations from the fit curve due to the 
failure of rotational invariance on a discrete lattice. In 
fact, fitting all the available data to Eq. ((8]) results in 
an increase of the reduced chi-squared without affecting 
appreciably the fit parameters. In Table I the results 
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cooling 


<t> 
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xi 


5 


0.886( 89) 


0.829(267) 


0.590(452) 


0.512(400) 


0.3 
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1.214( 66) 


0.782(145) 


0.567(270) 


0.485(365) 


0.2 
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1.590( 53) 


0.749( 82) 


0.538(162) 


0.453(325) 


0.2 
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1.998( 49) 


0.711( 51) 


0.530(110) 


0.444(314) 


0.1 


9 


2.423( 49) 


0.685( 39) 


0.508( 82) 


0.420(285) 


0.1 


10 


2.680( 23) 


0.652( 22) 


0.530( 44) 


0.444(313) 


0.3 



TABLE I. Summary of the fit values for SU(2) at (3 = 2.52. 

of our fit of Eq. jSJ to the SU(2) data at /3 = 2.52 are 
reported. The Ginzburg-Landau parameter k has been 
obtained through Eq. (|TU)) . 

In Figs. [2j [3l [4j [5]we display the fitted parameters versus 
the cooling steps. As regards the parameters /z, A/£„, 
and k, a short plateau is always visible. This corrobo- 
rates our expectation that the long range physics is un- 
affected by the cooling procedure. On the other hand, 
Fig. [2] shows that the overall normalization of the trans- 
verse distribution of the longitudinal chromoelectric field 
is more affected by the cooling. In fact, the parameter 4> 
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FIG. 2. SU(2): 4> versus cooling. 
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FIG. 3. SU(2): fj, versus cooling. Here and in other figures 
below data have been slightly shifted along the horizontal axis 
for the sake of readability. 

seems to displays an approximate short plateau after 9 - 
10 cooling steps in accordance with previous studies [l6j . 

B. SU(3) data 

We re-analyzed the SU(3) lattice data presented in 
Ref. [35j]. We have fitted the longitudinal chromoelec- 
tric field transverse distribution to Eq. © for (3 — 
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FIG. 4. SU(2): \/£ v versus cooling. 
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FIG. 6. SU(3): Ei at /3 = 6.0 after 10 cooling steps. 
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FIG. 5. SU(2): k versus cooling. 



5.9,6.0,6.05,6.1 and up to 16 cooling steps. Again, we 
find that Eq. ((8]) accounts for the transverse distribu- 
tion of the longitudinal chromoelectric field in the whole 
region x t > 0. In Fig. [5] we display our data for the 
transverse shape of the longitudinal chromoelectric field 
between static quark- ant iquark sources after 10 cooling 
steps at (3 = 6.0 together with the fit to Eq. ©. As for 
the SU(2) case we also display the points calculated at 
noninteger distances and checked that the fit to all the 
available data of Eq. ([SJ does not change the values of 
the fit parameters. 

In Table II we collect the results of our fit of Eq. © 
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4.564( 14) 


0.726( 8) 
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0.428(295) 
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0.449( 7) 
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7.644( 9) 


0.652( 4) 


0.412( 6) 


0.316(177) 


2.2 


8 


8.107( 8) 


0.630( 3) 


0.384( 5) 


0.287(150) 
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8.254( 8) 


0.612( 3) 


0.360( 4) 


0.263(130) 


5.5 


10 


8.227( 7) 


0.598( 3) 


0.341( 4) 


0.244(114) 
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8.108( 7) 


0.585( 3) 


0.325( 4) 


0.229(103) 


4.5 


12 


7.943( 7) 


0.574( 3) 


0.313( 4) 


0.217( 94) 


3.6 


13 


7.759( 6) 


0.563( 3) 


0.304( 4) 


0.208( 88) 


2.8 


14 


7.570( 6) 


0.553( 3) 


0.296( 4) 


0.201( 83) 


2.1 


15 


7.383( 6) 


0.544( 3) 


0.291( 4) 


0.196( 79) 


1.5 


16 


7.204( 6) 


0.534( 3) 


0.287( 4) 


0.192( 77) 


1.1 



TABLE II. Summary of the fit values for SU(3) at P = 6.0. 



to the SU(3) data at = 6.0 for cooling steps ranging 
from 5 up to 16. In Figs. M OH we show the fitted 
parameters versus the cooling steps. 
In fact, we see that the parameters /i, A/£„, and k, display 
a short plateau during the controlled cooling procedure 
as in the SU(2) case. Moreover, Fig. [7] shows that, at 
variance with the previous case, even the overall normal- 
ization of the transverse distribution of the longitudinal 
chromoelectric field <f> seems to displays an approximate 
plateau after 7-9 cooling steps in agreement with the re- 
sults of Ref. (ai. 



III. PENETRATION AND COHERENCE 
LENGTHS 

In Refs. [HI, Hi| it was found that the inverse penetra- 
tion length /i exhibits approximate scaling with the string 



5 




0.75 



2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 

cooling step 

FIG. 7. SU(3): <j) versus cooling. 
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FIG. 8. SU(3): n versus cooling. 



tension a. To check the scaling of our new determination 
of /i with the string tension, we use a parameterization for 
the SU(2) string tension obtained by means of a Cheby- 
shev polynomial interpolation to the string tension data 
collected in Table 10 of Ref. [E3|. 

In Fig.[TT]we display our determination of the ratio fi/y/a 
for three different values of j3. For comparison we also 
report /i/y/a, where the inverse of the penetration length 
\i is obtained by fitting the transverse profile of the lon- 
gitudinal chromoelectric field to Eq. (0| after 8 cooling 
steps (for details, see Ref. [Hj]). We see that our new de- 



0.5 



0.25 - 



"1 [T 



"i i r 



"i i r 



"i i r 



o 


(3=5.9 


□ 


p=6.0 


o 


P=6.05 


A 


p=6.1 



® 



J I L 



J I L 



J I L 



2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 

cooling step 

FIG. 9. SU(3): A/£„ versus cooling. 
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FIG. 10. SU(3): k versus cooling. 



termination of uJ_Ja is in satisfying agreement with the 
results of Ref. [35|. Thus, we confirm that /i displays an 
approximate scaling with the string tension a. Fitting 
our data for \ij yfa with a constant, we estimate 



\ij\fa — 4.133(98) 



(11) 



where the quoted error take care also of the systematic 
errors due to the scaling violations displayed by our data. 
Assuming y/a = 420 MeV, Eq. (fiT]) gives for the pene- 
tration length 
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FIG. 11. SU(2): n/y/5 versus /3. Open circles corresponds to 
the fit with Eq. Q after 8 cooling steps, full squares corre- 
spond to the fit with Eq. ([§} after 10 cooling steps. 




FIG. 12. SU(2): k versus /3 for 10 cooling steps. 



Moreover, we have also checked the scaling of 
the Ginzburg-Landau parameter k obtained through 
Eq. ([TO]) . In fact, Fig. [12] shows that k is almost in- 
sensitive to /3. By fitting the data with a constant we 
get 



0.467 ± 0.310 



(13) 



We would like to stress that our continuum extrapolation 
for the penetration length and the Ginzburg-Landau pa- 
rameter are in reasonable agreement with the results ob- 
tained in Refs. [H, EH HH . ^ n particular, we mav confirm 
that the Ginzburg-Landau parameter is consistent with 
the critical value k c = -4=, i.e. the SU(2) vacuum behaves 
as a dual superconductor which lies at the borderline be- 
tween the type I - type II superconductor regions. 



Also for the SU(3) gauge theory we studied the scaling 
of the "plateau" values of /i with the string tension. For 
this purpose, we have expressed our values of /i in units 
of i/ct, using the parameterization 



^G?) = /su ( 3)(ff 2 )[l + 0.2731 a 2 (g) 



(14) 



0.01545 a 4 (g) + 0.01975 a 6 ( 5 )]/0.01364 , 



a(g) 



fsv(3)(g 2 ) 



/su(3)(ff 2 (/3 = 6)) 



5.6 < @ < 6.5 , 



/su(3)(s 2 ) = (bog 2 ) 
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(15) 
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h = 
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given in Ref. [5^] . Figure [T3] suggests that the ratio [ij \fo 
displays a nice plateau in /3, as soon as (3 is larger than 
6. Accordingly, fitting the ratio fi/ ' \fu to a constant we 
get: 



2.799(38) , 



(16) 



which, assuming again the standard value for the string 
tension %fa = 420 MeV, corresponds to: 



A = — = 0.1676(23) fm . 



(17) 



The quoted error takes into account our estimation of sys- 
tematic effects due to the small scaling violations present 
in our lattice data. It is interesting to compare the 
present determination of the SU(3) penetration length 
with the one obtained previously [35j by fitting the lat- 
tice data to Eq. (QJ. In fact in Ref. [35[ we obtained 
^= = 2.325(5) which, unlike to the SU(2) gauge theory, 

seems not to agree with Eq. (jT6j) ■ This discrepancy must 
be ascribed to the fact that, as discussed below, in the 
present case the coherence length exceeds the penetration 
length. In other words, the SU(3) vacuum behaves like a 
type I superconductor and, as we already discussed, the 
London equation Eq. Q is not good enough to describe 
the transverse distribution of the longitudinal chromo- 
electric field. 

As concerns the Ginzburg-Landau parameter, in Fig. [14] 
we present our lattice data for different values of j3. In- 
deed, k is almost insensitive to j3. By fitting the data 
with a constant we get 



k = 0.243 ± 0.088 



(18) 



which confirms that k < k c (type I superconductor). It 
is worthwhile to note that our Eqs. ([T7]) . (|18p are not in 
agreement with the recent determinations in Ref. [53| . 
where it is reported A = 0.2013(174) fm and k = 
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1.218(109). We believe that the origin of the discrep- 
ancies resides in the use of different lattice operators to 
extract the longitudinal chromoelectric field. In fact, the 
authors of Ref. [53] use a lattice operator which is sen- 
sitive to the square of the chromoelectric field instead 
of our correlation function, Eq. ([TJ , which measures the 
chromoelectric field strength. In any case, we believe 
that these discrepancies deserve further studies. To 




6.05 6.1 6.15 



FIG. 13. SU(3): n/yfa versus /3. Full squares correspond to 
the fit with Eq. ((8} after 10 cooling steps. 
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different study of SU(2) and SU(3) vacuum in a constant 
external chromomagnetic background field 541 ] . In fact, 
in Ref. [HH numerical evidence was presented that the 
deconfinement temperature for SU(2) and SU(3) gauge 
systems in a constant Abelian chromomagnetic field de- 
creases when the strength of the applied field increases. 
Moreover, as discussed in Refs. [Ill, [54, 55], above a criti- 
cal strength \JgH c of the chromomagnetic external back- 
ground field the deconfined phase extends to very low 
temperatures. It was found [54j | that the ratio between 
the critical field strengths for SU(2) and SU(3) gauge 
theories was 



SU(2) 



Vff^c|sU(3) 



2.03(17) 



(19) 



It is interesting to compare the ratio between the critical 
field strengths Eq. (fl9| with the analogous ratio between 
penetration and coherence lengths. Combining Eqs. ([S]), 
([IT ]) . ([15 ]) . (TT6"]) . and tHHJ), we readily obtain 
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FIG. 14. SU(3): k versus /3 for 10 cooling steps. 



It is remarkable that the ratio between the penetration 
lengths, respectively for SU(3) and SU(2) gauge theories, 
agrees with the analogous ratio between the coherence 
lengths, albeit within the rather large statistical uncer- 
tainty. Moreover, both ratios are in fair agreement with 
the ratio between the critical field strengths, Eq. ([19]) . 
As stressed in the Conclusions of Ref. [54] . the pecu- 
liar dependence of the deconfinement temperature on the 
strength of the Abelian chromomagnetic field gH could 
be naturally explained if the vacuum behaved as a disor- 
dered chromomagnetic condensate which confines color 
charges due both to the presence of a mass gap and the 
absence of color long range order, such as in the Feynman 
picture for Yang-Mills theory in (2+1) dimensions [56j . 
The circumstance that the ratio between the SU(2) and 
SU(3) penetration and coherence lengths agrees within 
errors with the above discussed ratio of the critical chro- 
momagnetic fields, suggests us that the Feynman picture 
of the Yang-Mills vacuum could be a useful guide to un- 
derstand the dynamics of color confinement. 



IV. CHROMOELECTRIC FLUX TUBE AND 
STRING TENSION 



summarize, we have found that the transverse behavior of 
the longitudinal chromoelectric field can be fitted accord- 
ing to Eq. @ for both SU(2) and SU(3) gauge theories. 
This allows us determine the coherence and penetration 
lengths. In Ref. 35] it was stressed that the ratio between 
the penetration lengths respectively for SU(2) and SU(3) 
gauge theories recalls the analogous behavior seen in a 



We have shown 16, 35] that the color fields of a static 
quark-antiquark pair are almost completely described 
by the longitudinal chromoelectric field, which in turn 
is approximately constant along the flux tube. This 
means that the long-distance potential acting on the color 
charges is linear. Using our data and the parameteriza- 
tion Eq. for the chromoelectric flux tube, we are able 
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FIG. 15. SU(2): y/a/ii versus /3. Full points correspond to 
Eq. (|24|l ; open points refer to the lattice string tension (see 
the discussion in the text). 
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FIG. 16. SU(3): \/<j/fJ. versus /3. Full points correspond to 
Eq. (|24|l ; open points refer to the lattice string tension (see 
the discussion in the text). 



to compute the string tension given as the energy stored 
into the flux tube per unit length: 



cr Bi ~ i / ,l-.r l E'f(.r l ) . 



(22) 



where, to avoid confusion, we have denoted the flux-tube 
string tension as er^,, while a will indicate the lattice 
string tension. It is worth to note that the string ten- 
sion (jEi defined by Eq. (j2"2")l does not depend on xi as 



long as the longitudinal chromoelectric field is constant 
along the flux tube. Obviously this last condition is not 
strictly fulfilled on a finite lattice. From Eqs. (|22j) and (jHJ) 
we obtain an explicit relation between the string tension 
and the parameters 4>, fi, and a of the fit Eq. ([8]) to the 
chromoelectric flux tube profile: 



2 ,,4 



<JE, 



1 0V 



1 



4tt a 2 Kf(a) Jn 
After performing the integration, we get 



drrK 2 ((^ 2 r 2 + a 2 ) 1 ' 2 ) 



8tt 



1 - 



KM 
K 2 (a) 



1/2 



(23) 



(24) 



Naively, one expects that the string tension defined as 
the energy per unit length stored into the flux tube chro- 
moelectric field Eq. (I2~il) should agree, at least approx- 
imatively, with the string tension measured on the lat- 
tice. To check this, in Figs. (ITS"]) and ([T6|) we compare 
the flux-tube string tension in Eq. (|24p with the lattice 
string tension (obtained as detailed in the previous Sec- 
tion) for SU(2) and SU(3) gauge theories, respectively. 
It is evident from Figs. (|T5|) and (fT6|) that 



(25) 



for both SU(2) and SU(3) gauge theories. At first sight 
this result looks quite surprising. In fact, the lattice 
string tension should contain the total energy per unit 
length stored into the flux tube. As a consequence we 
can write 



cte, 



^cond 



(26) 



where er conc i takes into account the contribution due to 
the order parameter condensate. We may, in turn, ob- 
tain an estimate of this contribution to the total string 
tension as follows. Since within the vortex core the order 
parameter condensate vanishes, we have 



^cond 



?r£ £ C ond 



(27) 



where e con d is the condensation energy density and £ is 
approximately the vortex core size. Note that the minus 
sign is due to the loss of condensation energy in the nor- 
mal region with respect to the confining vacuum where 
the order parameter is nonzero. Now, it is usually as- 
sumed that in the confining vacuum it is energetically 
favored to have a condensation of the order parameter as 
in ordinary BCS superconductors where the supercon- 
ducting transition is energetically driven by the coher- 
ent condensation of Cooper pairs. Therefore, one is led 
naturally to suppose that e con d < 0. Thus we see from 
Eqs. (f26| and ((27|) that one should obtain cte, < a. On 
the contrary our numerical results clearly indicate that 
a Ei exceeds a. The only possible conclusion we can de- 
rive is that the order parameter condensation energy is 
positive. Thus, the confining transition must be driven 
by disordering the gauge system. In other words, even 
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though the condensation of the confining order param- 
eter costs energy there is a huge number of degenerate 
physical configurations such that the configurational en- 
tropy easily overcomes the energy cost. This means that 
the deconfining transition is an order-disorder transition, 
much like the Bcrczinskii-Kostcrlitz-Thouless transition 
than the BCS superconducting transition. It is remark- 
able that this conclusion reinforces our previous picture 
of the confining vacuum that behaves like a disordered 
chromomagnetic condensate which confines color charges 
due both to the presence of a mass gap and the absence 
of color long r ang e order, such as in the Feynman quali- 
tative picture [561 ] . 

V. CONCLUSIONS 

In the present paper we studied the chromoclectric 
field distribution between a static quark- antiquark pair 
in SU(2) and SU(3) pure gauge theories. By means of 
the connected correlator given in Eq. (fTJ) we were able to 
compute the chromoelectric field that fills the flux tube 
along the line joining a quark-antiquark pair. Since our 
connected correlator is sensitive to the field strengths in- 
stead of the squared field strength, we were able to fol- 
low the transverse shape of the color fields up to sizable 



distances. Using some dated results in ordinary super- 
conductivity based on a simple variational model for the 
magnitude of the normalized order parameter of an iso- 
lated vortex, we proposed that the transverse behavior 
of the longitudinal chromoelectric field can be fitted ac- 
cording to Eq. ([5]). which allowed us to get informations 
on the penetration and coherence lengths. In fact we 
found that our Eq. (JS) is able to reproduce the trans- 
verse distribution of the longitudinal chromoelectric field 
in the whole available region. In the case of the SU(2) 
gauge theory we argued that the confining vacuum be- 
haves as a dual superconductor which lies at the border- 
line between the superconductor type I - type II regions. 
On the other hand, we found that the SU(3) vacuum be- 
longs to the superconductor type I region. We found that 
the ratio between the penetration lengths respectively for 
SU(3) and SU(2) gauge theories agrees with the analo- 
gous ratio between the coherence lengths, albeit within 
the rather large statistical uncertainty, and both ratios 
are in fair agreement with the ratio between the criti- 
cal chromomagnetic fields. Finally, we suggested that 
the deconfining transition resembles the order-disorder 
Berezinskii-Kosterlitz-Thouless transition and that the 
confining vacuum behaves like a disordered chromomag- 
netic condensate in agreement with the Feynman quali- 
tative picture of the Yang-Mills vacuum. 
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